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SUMS AND DIFFERENCES OF POWER-FREE NUMBERS 


JULIA BRANDES 

Abstract. We employ a generalised version of Heath-Brown’s square sieve in 
order to establish an asymptotic estimate of the number of solutions a, b to the 
equations a + b = n and a — b = n, where a is fc-free and b is 1-free. This is the 
first time that this problem has been studied with distinct powers k and 1. 


1. Introduction 

Let Afc C U denote the set of fc-free numbers. We are interested in the number 
R^i{x] n) of twins a,a — n<x with a given distance n, where a is k-iree and a — n 
is /-free. A closely related question asks for the number R^i{x]n) of solutions 
a & Ak and /> G A; to the equation a + b = n, where we may impose the additional 
size restriction that one of the summands be bounded by a parameter x. Observe 
that trivially R^i{x]n) = 0 when n> x and Rli{x;n) = R'^i{n;n) when x > n. 
More precisely, a is restricted to an interval of length X^{x] n), where X~{x] n) = 
max{0, x — n\ and A+(a;; n) = min{a;, n}, so as the k-iree numbers are of positive 
density in the integers, we expect the main term to be proportional to X^{x]n). 

When the exponents k and I are equal these are well-studied problems, and it 
is not hard to establish an asymptotic formula. In fact, if u{k) denotes the least 
number u satisfying 

n n (1-2 pT +O ((! + «)"+'), 

p prime ^ p prime 

p^\n 

then elementary methods (see Carlitz [3] for R 2 2 {x] 1), and Evelyn and Linfoot [6] 
and Estermann [5] for R'l^^{n]n)) suffice to show uj{k) <2/{k + 1). Interestingly 
enough, this still seems to be the state of the art in the case of while 

for R'^j^iyX] 1) better error terms are available. In fact, Heath-Brown |7j used sieve 
methods to establish oj{2) < 7/11, an approach that has later been generalised to 
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higher powers in the author’s diploma thesis jT], yielding an error terirQ of 

if A: < 4, 


oj{k) < 


14/(7A; + 8) 
8A; + 3 


if /c > 5. 


+ 6A; + 2 

The advent of the determinant method triggered further improvements. Dietmann 
and Marmon |1] established u:{k) < 14/(9fc), while independently Reuss [10] re¬ 
duced the exponent to u{2) < (26 -|- \/433)/81 0.578 and u{k) < 169/(144fc) 

for larger values of k. Both results are much stronger than [Ij. 


It is now natural to ask whether these methods continue to be applicable when 
the powers are distinct. However, except for a passing remark by Briidern et 
ah [2] claiming that their approach to study via the circle method 

can be adapted to distinct powers, the author is not aware of any previous work 
regarding this more general setting. In this note, we show how the power sieve as 
developed by Heath-Brown [7] in the case k = 2 and extended to higher powers 
independently by the author [T| and Munshi [9] can be adapted to obtain an 
asymptotic formula for Rj^i{x;n). Let 

n ( /-7'0 “d 

p prime ^ ^ / p prime 

p'^ln 

then the result is the following. 


Theorem 1. Suppose 2 < k < I, then the counting functions Rffx^n) and 
R'^i{x\n) obey the asymptotic formula 

Rfi{x;n) = X^{x;n)Ck,i{n)&k,i + 0{{x + 


where the exponent u{k, 1) in the error term is bounded by 


A{k + l)+3 


u{k, 1) < 


2{2k -|- 1)(/ -|- 1) 
7{k + l) 

Qk + 2l + 7kl 
k + l 


, k(^l -\- 2 ) 


if k <l < {A/3){k — 1), 
if {A/3){k — !)</< 4A;, 


if I > 4fc, 


and the implied constant depends on k, I and e, but is independent of x and n. 
In the case k > I the same asymptotic formula is true with the roles of k and I 
reversed. 


^The original work contains an oversight in the numerical analysis, which led to a different 
exponent. This error is corrected in the present paper. 
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Observe that Theorem [U reproduces the result from [T] and extends it to the 
problem of estimating n), in which case it is the hrst improvement on the 

error term since Estermann’s work from 1931 |5]. For x n the error remains 
non-trivial as long as -c x. 

A word is in order as to why we believe this result to be competitive, in par¬ 
ticular in the light of the strong bounds that can be obtained by the determinant 
method in the case k = 1. Observe that for I > k every fc-free number is a forteri- 
ori also /-free, so when / is not much larger than k the bounds on u:{k) obtained by 
Reuss ra will yield a stronger error term than that of Theorem [T] However, in its 
current formulation the determinant method fails to properly accommodate dis¬ 
tinct powers as it seems to rely crucially on rewriting the equation u^h — v^j = n 
in the shape {u/vY — j/h = n/{v^h) and then counting rational points within a 
distance On{{v^h)~^) from the curve a transformation impossible when 

k and I are distinct. This means that when I is sufficiently large in terms of k 
we should expect Theorem [1] to yield a stronger error term than what is being 
obtained by the determinant method, and indeed, our exponent supersedes that 
of Reuss dni as soon as I is larger than roughly 144fc/25. Furthermore, in the limit 
I —)■ oo the exponent in Theorem [T] tends to 1/k, thus replicating the known error 
term in the asymptotic formula for the number of /c-free integers. For comparison, 
the methods of Briidern et al. j2] yield an error term whose exponent is bounded 
below by 1/2. 

Throughout the argument the following conventions will be obeyed. No effort 
will be made to track the exact ‘value’ of inhnitesimal quantities arising in the 
analysis, so all statements involving e are true for any e > 0. The implied con¬ 
stants will be allowed to depend on all parameters except n and x. We write 
d{n) for the number of divisors of an integer n, and (a, h) denotes the greatest 
common divisor of two integers a and b. Furthermore, since the analysis differs 
only very little for and R~ we will write and use the ambiguous signs ± 
and =F throughout. This will happen in a consistent way, so in =p the negative 
sign corresponds to R'^ and the positive value arises when R~ is considered. 

The author is grateful to the referee for valuable comments. 


2. Elementary considerations 
Recall that our counting function is given by 

R^i{x] n) = {a < X ■. a G Ak, ±(n — a) G Ai}. 
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Using the indicator function on /c-free numbers jJikin) = 
rewrite this as 


^^{u)^i{v)N^^{x,n■u,v), 

u,v<x 


where 


N^i{x, n; u, v) = Card{l < a < x : ±(n — a) > 1, u^\a, v^\{n — a)}. 

One sees that with our assumption k < I one has N^i{x, n; m, n) = 0 when (m, 
does not divide n and 

Nk,i{x,n-u,v) = + Q(l) 

otherwise. 

For a parameter y to be chosen optimally later the contribution with uv <y is 
given by 


X^{x-n) 

w,i;EN 
{u^v)^\n 


IJ,{u)fi{v){u, vY 

ykyl 



+ o(Ei)^ ( 1 ) 

uv<y 


We may rewrite the sum in the main term by extracting common factors and 
writing u = u'd and v = v'd. Note that fi{u'd) takes a non-zero value only if u' 
and d are coprime, so we have 

E jj,{u)y{v){u,vY y{u')y,{v') 

y^kyi 2^ yi 2-^ u'^v'^ 

{u^v)^\n d^\n {u'^v')=l 

(d^u' )={d^v' )=\ 

d} 

df‘\n {h,d)=l w\h 

The sum over h is multiplicative with an Euler product expansion given by 

n(i - - p"')= &k,i n(i - p"^ - p"')”\ 

p\d p\d 

and after expanding the d-summation one recovers the constant k < 1. 

The respective result for I > k follows by symmetry. 


We now turn to the error 


= RkMi^) - X^{x]n)&k,iCkM)- 
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Using that X^{x]n) < x, we have 

^ 1+x ^ Ku)K'v)N^,i{x,n-u,v), 

Uj u 

uv<y uv>y u,v<x 

{u,v)'^\n uv>y 

where the first two terms are the error terms from ([T]) and the last term acconnts 
for the contribution not considered in ([T]). Observe that the first term is 0{y\ogy) 
trivially and the second term is bounded above by 

uv>y d>^\n v=l mSN 

{u,v)^\n u>yl{(Pv) 


The inner sum in the above expression is {(Pv/yY ^ if n < y/(P and -C 1 else, 
whence on the right hand side we have 


fj^k\n v=l 1/GN 

u>y/{d?v) 


u ^ xy 


d^\n 


k-l-l 


d^ln «eN 
v>yld? 


-I 


The n-summation in the first term on the right hand side is empty unless d < 
and in the second term we have another case distinction as before according to 
whether y/d?' > 1 or not. It follows that the above is 


•c 

xy^~ 

-f xy^~ 


+ X Y^ d. 



d<Vy 

d<Vy 

d>^ 



d'^\n 

dl^ln 

d}^\n 

< 

xy^~ 


xy~^^'^rf 


< 

xy^~ 





This yields the overall error 

E^i{x] n) -C y^^^^ + xy^~’^^^ + xy~^l‘^rf + ^ y{u)y{v)N^i{x, n; u, v). 

uv>y 

u,v<x 

The last term will be the focus of our analysis in the following section. 


3. Sieving for powers 

We split the terms with uv > y in roughly (logx)^ intervals U < u < 2U and 
V < V < 2V. Then there exist values U and V such that 

f^Nf^(v)X^i(x,n;u,v) < W(loga;)^ 

uv>y 
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where N is given by 

N = Card{M, v,h,j: U < u < 2U, V < v < 2V, hv!^ ± jv^ = n}. (2) 

Let Nh denote the number of n, v,j counted by N for a hxed value of h, and write 

^ N,. 

H<h<2H 

There are (9(logx) intervals H < h < 2H with 

H < xU-\ (3) 

so for some H one has N -C N{H) logo: and consequently 

Efi{x]n) <C + xy~’‘^‘^rf + (loga:)^iV(iL). (4) 

We may assume without loss of generality that 

U^>VK ( 5 ) 

This might involve swapping the parameters h, u and k with j, v and I, respec¬ 
tively, which in the situation of counting a — b = n inserts a minus sign, so we let 
(T = — 1 if this sign change occurs and a = +1 else. Then, having hxed a suitable 
value H, let 

Ji = min{l, ±2-^V-^{an - 2^+^i7[/^)}, J 2 = ±V'(cTn - 
Note that by construction we have Ji < J 2 . It follows that 

Nh < Card{ Ji < j < J 2 -, V <v< 2V, U < u < 2U : hu^ ± jv^ = an}. 


Let Wk{z) = 0 for all 2 ; that are not divisible by and 
Wk{h^~^a) = Card{(n, j): h\{an =F v'‘j), a = cm =F 


where V < v < 2V and < j < ^ 2 - This is legitimate as an =|= v^j > 0 for all 
admissible values, and we have 


Nh < 

m=l 


We now apply the power sieve as in jl]. Let Q be a parameter satisfying logx < Q < x 
which will be chosen optimally later, and let V be the set of primes in the inter¬ 
val [Q, 2Q] that do not divide hn and are congruent to 1 modulo k, then by the 
Siegel-Walhsz Theorem we have Card V zz Q/ log Q. Write further for a non¬ 
principal character modulo p whose fc-th power is the principal character modulo 
p. That such a character exists is ensured by our choice of the set V. The power 
sieve as in [9l Lemma 2.1] or [H Theorem 1] now yields 


Nh<^ 


logQ 

Q 


Z=1 


wJz) + max 
p,q&'P 
p¥=g 






^{an ^ v’'j)\ f h’‘~^{an ^ v’-j) 


p 


, ( 6 ) 
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where the sum is over aXW < v < 2V, < J < J 2 satisfying h\{an =F v^j). The 

hrst term can be estimated by 


logQ 

Q 


cx> 

'^Wk{z) < 

z=l 


^ V<v<2V ^<x 

^=<jn (mod v^) 


so the deeper challenge lies in establishing a good bound for the second term. 
This analysis will follow the treatment of [7] and jT] closely. 

Denote the expression inside the absolute values in the second term of (jSD by 
S. We have 


hpq 


S = (hpq) ^ ^ S{h,pq]'y,6)e-y(j)s, 

7,(5=1 


where 


hpq 


S{Kpq-,'y,S) = 


a,0=1 
h\{an^a^ 0) 


an =F a^/3 \ / an =F a'-jS \ f 'ya + 6(3 

e 


p 


q 


hpq 


E 


V<v<2V 


E - 


—yn 
hpq 

-5j 

hpq 


min V, 


7 


hpq 


-1' 


< min V-\U'^H + n), 


hpq 


-V 


Writing further 


„ / f crn^ a’-(3 \ f ca + df3 

Si{p-,c,d) = 2 _^ I - I e 


a,13=1 


S 2 {r‘-,c,d)= E 


P /t 

ca + djS 


p 


Q:,/3=1 
r'^\{an^^o0 0) 


(7) 

( 8 ) 

(9) 


one has 


S{h,pq;'y,5) = Si{p; c, d)Si{q; -c, -d) JJ S' 2 (r*;c,d) 

r'^\\h 

where h = Hr* denotes the prime power decomposition of h, and c and d are 
certain integers satisfying (c, hpq) = ( 7 , hpq) and (d, hpq) = {6, hpq). 
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Lemma 1. Suppose p is a prime number not dividing n, then the exponential 
sums Si{p;c,d) and S 2 {p;c,d) can be estimated by 

Si{p;c,d) -Cp, 

S 2 {p; c, d) -C p^^'^{p, c, 

Proof. This is analogous to the treatment in [71 §5] or jU Lemma 2] and is ele¬ 
mentary except in the generic case when p \ cd, in which case it depends on the 
Weil conjectures. □ 

Let w denote the product of all primes coprime to n that divide h exactly once. 
By following through the arguments of [71 p.256sq.] or [H p.l2] one sees that 
Lemma [T] implies 

JJ S' 2 (r*; c, d) -C 7 , 

r*||/i 

whence we obtain 

S{h,pq;^,6) < 

Inserting the definitions (jT]), ([HD and ([HD of S, 9.y and (fs, respectively, together 
with dHD, one concludes that 

S <^Q-^H-^w-^/^+^(v^-\x + n)w^/^+ VQ^H ^ 

\ l<S<hpq/2 

+ V-^Q^H{x + n) 7~^(M;,7)^/2 + Q^Lr2^(75)"^(u;,7,(5)^/M 

l<7<hpg/2 7,5 / 

< + n)w^ + {V-\x + n) + V + HQ^) uj-^/^+%\og HQ^. 

We can now estimate our sifted sum N{H). By an argument similar to that in 
m p.257] or [H p.l3] one has YliH<h<2H'^ gQ is bounded 

by 

(a; + n)Q-V^-' + {y-\x P n) PV P HQ^)EY- 

After recalling ([HD, it follows that 

N < Q-W^-\x P n)^+" + n)U-'^/^V-^ 

^ ^3/2+£fj-3fc/2g2_ gg) 


It remains to optimise the parameters in order to establish the theorem. 
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4. A GAME OF OPTIMISATION 


In this section it is our task to optimise the parameters in the error 

n) <C + xy~'‘^‘^rf + x^N, 

where N is bounded by the expression in (ITOl) . We will write z = x + n and use 
X z throughout. With this convention the hrst and the last term of (ITOl) are 
roughly equal if 

Q = (loga:)2 + ^-i/6+£f/fc/2p.(i-o/3_ 

Furthermore, it follows from the assumption made in ([3]) that 

ljjkYl/2y ^ ijj\ky{l-\)lyl/2y 

for every A G [0,1]. With A = {I — 2)/{k + l) the powers in U and V coincide, and 
since we had UV 3> y this expression is bounded by ^ yielding 

N -C ^F6+£f/-fc/2^2(l-/)/3 Yl2+ey-kl2y-l _j_ ^1/2+e^ . (H) 

In order to evaluate the terms depending on U and V one has to make recourse 
to an auxiliary bound for the quantity given in (j2|). 

Lemma 2. We have 


Proof. Note that by the dehnition ([2]) the quantity N is equal to 

V<v<2V h<xU-^ U<u<2U 

u^h=an (mod v^) 


( 12 ) 


The innermost sum can be estimated via Hensel’s Lemma. If d = (n,n^), the 
quantity in question is the number of u G {U, 2U] satisfying dlu’^ for which u^h/d = 
an/d (mod /d). This congruence has at most k solutions for each prime divisor 
of /d^ giving altogether an upper bound of dk{v). On the other hand, one has 
the trivial bound UV~K Inserting this into fll2p yields the required bound. □ 


Notice that one has 

xl+eul-kyl-l ^ ^ ^ 

which is absorbed in the second term of (jl]). We obtain our bound by comparing 
the remaining term of f[T2D with the hrst two terms in (ITT]) , respectively. Observe 
that 

min zU~^V^ ^ ^(3/2)a-|-(l-a)^-(fc/2)a-(l-a)fcpr-Zo-|-(l-a) 

= fcpA—(Z-l-l)a-l-l 


As before, it follows from ([5]) that for every y G [0,1] one has 

jjka/2—ky—{l+l)a+l ^ jj^ik{a/2—l)y{l—^i)l{a/2—l)—{l+l)a+l 
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and on choosing /i = powers in U and V coincide and one obtains 

Similarly, we have 

min |//6^-fc/2v^2(i-z)/3^ zU-^V} 

^ ^l+/3/6^fc^/2-A:-^-(2Z+l)/3/3+l 
< ^l+/3/6f/;^fc(/9/2-l) pA(l-i^)Z(/3/2-l)-(2/+l)/3/3+l 


and the powers of U and V coincide on choosing u = • The bound for 

N is therefore given by 


jY 2(^+0^ -)- ^^+“/2 +s^M^(“/ 2—1) _j_ y>^k{fi/2—1) 


(13) 


Put y = z^ for some suitable power a; < 1. This is legitimate whenever < x, 
and in the opposite case the statement of the theorem is trivial. Then by inserting 
the values for y and v we see that the last two terms of (1T3|) coincide if /9 = 3a, 
except possibly in the case when uj = ujq = {k + l)/{kl + 2k). However, one checks 
that in this case the two terms coincide for all choices of a and /5. Also, the last 
term in (IT^ is roughly of size y if 


6a;/(l + fc) — 6(A; +/) 

ujk{l + 2) — {k + 1) 

It remains to hnd the minimal u for which the parameters a,/3,/i, u all he in 
the unit interval [0,1]. Notice that our choice of (3 implies 

(2 + /)u;-3 

^~~{2k + l + 2kl)u-2{k + l)' 

{2 + l)uj-3 

{2k + 3Z + 4ikl)uj — 4(A; + 1) 


The horizontal asymptotes of /3, y and u as functions of oj are negative, so 
the pre-images of [0,1] under each of these functions are also closed intervals. A 
modicum of computation shows that the endpoints of these intervals are given by 


k + l 


— 

kl + V 

i^y=i = 

^11=0 = 

3 

= 

/-t-2’ 

^iy=0 = 

3 

(^1^=1 = 

/-f2’ 


7{k + l) 

2k + Ql + 7kV 
2{k+ !) + ?> 
2{k + l){l + iy 
4(/c 3 

2(fc -|- 1) (2/ -|- 1) 


In order to analyse the intervals dehned by these values, one has to consider 
various cases. The expressions for /3, y and v are increasing functions in u when 














SUMS AND DIFFERENCES OF POWER-FREE NUMBERS 


11 


I > 2k and decreasing else. Consider first the case when I > 2k. Then one has 
01^=1 > so the optimal valne for u is at the lower endpoint of the interval 

[oJv=Q.,OJy=]\ n [01/3=0, Ci;/3=i]. 

One checks that this is non-empty and has its lower endpoint at 01 / 3=0 = If, 
on the other hand, I is smaller than 2k, then one has 01^=1 < oij,=i and we want 
to hnd the lower endpoint of the interval 

[oij,=i, 01,y=o] n [oi/3=i, 01/3=o]. 


This, too, is always non-empty, and one has 00^=1 < ^^/ 3 =i if ^ < (3/4) A; -|- 1 and 
the opposite ineqnality whenever (3/4)fc -t-1 < / < 2k. Finally, in the case I = 2k 
all relevant valnes of 01 coincide and conseqnently there is nothing to optimise. 
Hence the error is given by with 


oii(A:,/) 


r 7{k + l) 

2k+ 61 + 7kl 
4(A; -|- Z) 4- 3 
2(^k + 1 ) ( 2 / -|- 1 ) 
k + l 
^ kl +1 


if I < {3/A)k + 1, 
if (3/4)A; + 1 < Z < 2fc, 
if Z > 2k. 


Recall now that dnring the analysis we possibly swapped the parameters k 
and Z. Reversing that swap, we may retnrn to the original assnmption that k < I 
and conclnde that the optimal exponent is given by 


UJ2{k,l) 


4(Zc -|- Z) -|- 3 
2{2k 4- 1)(Z 4- 1) 
7{k + l) 

6k+ 21 + 7kl 


iik<l< (4/3)(A;- 1), 
if (4/3)(fc- 1) < Z. 


Analogonsly, the hrst term of (fT3j) is optimised by ^i/ 2 +e^-hfc- 2 )/( 2 fc-i- 2 ) which co¬ 
incides with for y = ^(^+0/(2fc+fc0^ exponent exceeds U 2 {k,l) when¬ 

ever Z > 4Z^. The proof is now complete on noting that the remaining terms 
and ^i+S|/-V 2 jj/^ (gjj coincide with y when y = y = ^2/(i-i-2)^ respectively, 

both of which are dominated by the errors previonsly fonnd. 
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